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Abstract: 

We snggest an anomalons effective Lagrangian which reprodnces the anomalons conformal 
and chiral Ward identities and topological charge qnantization in QCD. It is shown that 
the large Nc Di Vecchia-Veneziano-Witten effective chiral Lagrangian is locally recovered 
from onr resnlts, along with 1/Nc corrections, after integrating ont the heavy “glneball” 
helds. All dimensionfnl parameters in onr scheme are hxed in terms of the qnark and 
glnon condensates and qnark masses. We argne that for a certain range of parameters, 
metastable vacna appear which are separated from the trne vacnnm of lowest energy by 
domain walls. The snrface tension of the wall is estimated, and the dynamics of the wall is 
discnssed. The U(l) problem and the physics of the psendo-goldstone bosons at different 
6 angles are addressed within the effective Lagrangian approach. Implications for axion 
physics, heavy ion collisions and the development of the early Universe dnring the QCD 
epoch are discnssed. 



1 Introduction 


The effective Lagrangian techniques have proved to be a powerful tool in quantum field 
theory. Generally, there exist two different definitions of an effective Lagrangian. One 
of them is the Wilsonian effective Lagrangian describing the low energy dynamics of 
the lightest particles in the theory. In QCD, this is implemented by effective chiral 
Lagrangians (ECL’s) for the pseudoscalar mesons, which are essentially constrained by 
the global non-anomalous SU{Nf) x SU{Nf) and (for large W) anomalous U{1) chiral 
symmetries. Another type of effective Lagrangian (action) is defined as the Legendre 
transform of the generating functional for connected Green functions. This object is 
relevant for addressing the vacuum properties of the theory in terms of vacuum expectation 
values (VEV’s) of composite operators, as they should minimize the effective action. 
Such an approach is suitable for the study of the dependence of the QGD vacuum on 
external parameters, such as the light quark masses or the vacuum angle 6. The lowest 
dimensional condensates (TT), (G^), (GG), which are the most essential for the QGD 
vacuum structure, are related to the anomalously and explicitly broken conformal and 
chiral symmetries of QGD. Thus, one can study the vacuum of QGD with an effective 
Lagrangian realizing at the tree level anomalous conformal and chiral Ward identities of 
the theory. The utility of such an approach to gauge theories has been recognized long ago 
for supersymmetric (SUSY) models, where anomalous effective Lagrangians were found 
for both the pure gauge case and super-QGD (SQGD) p. 

The purpose of this paper is a detailed analysis of the anomalous effective Lagrangian 
for QGD with Nf light flavors and W colors (more precisely, of its potential part) which 
was suggested recently by two of us, and briefly described in |^ . It was obtained as a gen¬ 
eralization of an anomalous effective Lagrangian for pure YM theory which was proposed 
earlier in [Q. The constructions of ^ and @ can be viewed as non-supersymmetric coun¬ 
terparts of the Veneziano-Yankielowicz (VY) effective potential [Q] for SUSY YM theory 
and the Taylor-Veneziano-Yankielowicz (TVY) effective potential for SQGD, respectively. 
The results obtained in ^ reveal some striking similarities between the supersymmetric 
and non-supersymmetric effective potentials and the physics that follows. Notably, the 
effective potentials for QGD and gluodynamics are holomorphic functions of their fields, 
analogously to the SUSY case[|. Moreover, they have both “dynamical” and “topological” 
parts - a structure which is similar to that of the (amended |^) VY effective potential 
1]^. As will be discussed in detail below, this “topological” part of the effective potential 
turns out to be crucial for the analysis of the physical 6 dependence in QGD. 

The interest in such an effective Lagrangian for anomalously broken conformal and 
chiral U{1) symmetries is several-fold. First, it provides a generalization of the large 
Di Vecchia-Veneziano-Witten (VVW) EGL ^ (see also [0) for the case of arbitrary Nc 
after integrating out the massive “glueball” fields. At this stage, no holomorphy is present 
in the resulting effective chiral potential. In this way we arrive at a Wilsonian effective 

^It should be noted that the status of holomorphy in SUSY and ordinary QCD is different. In 
supersymmetric theories, holomorphy is a consequence of supersymmetry and, moreover, holomorphic 
combinations zLiGG are determined by the structure of the anomaly supermultiplet. On the contrary, 
holomorphy in the ordinary YM theory (in a special sense) follows provided we make some plausible 
assumptions which are shown to be self-consistent a posteriori, see Sect.2. 
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Lagrangian for the light degrees of freedom consistent with the Ward identities of the 
theory and a bnilt-in qnantization of the topological charge (see below). One may note 
that in principle snch a Wilsonian effective Lagrangian satisfying all general reqnirements 
on the theory could also be written down without constructing hrst a more complicated 
anomalous effective potential including also the “glueball” degrees of freedom. The result¬ 
ing effective potential is found to contain a sine-Gordon term whose large W expansion 
reproduces the VWW effective potential in the vicinity of the global minimum, along with 
1/W corrections. The presence of such term in the effective potential implies that the 
theory sustains the domain wall excitations. This observation may be important in the 
contexts of cosmology and heavy ion collisions. Furthermore, our two-step approach to 
the derivation of the effective chiral Lagrangian has an additional merit in that all dimen¬ 
sionful parameters in our scheme are hxed in terms of the gluon and quark condensates 
and quark masses. (The only entries which are not hxed in our scheme are two dimen¬ 
sionless integer-valued parameters related to the vacuum structure and 9 dependence of 
the theory. As their values are still a subject of some controversy, see Sect. 3, in most 
of the paper we will keep them as free parameters.) The absence of free dimensionful pa¬ 
rameters helps to better understand the origin of the rj' mass (the famous U(l) problem). 
In particular, it yields a new mass formula for the rj' for hnite W in terms of quark and 
gluon condensates in QCD (see Eq.(^) below). Second, it allows one to address related 
questions of the phenomenology of pseudoscalar mesons, such as vr^ — 77 — 77 ' or 77 — 77 ' mixing 
with no further phenomenological input. Third, such an effective Lagrangian allows one 
to address the problem of 6 dependence in QCD. In contrast to the approach of Ref.|^ 
which deals from the very beginning with the light chiral degrees of freedom and explicitly 
incorporates the 17 a ( 1 ) anomaly without restriction of the topological charge to integer 
values, in our method both the 17a(1) anomaly and topological charge quantization are 
included in the effective Lagrangian framework. After the “glueball” helds are integrated 
out, the topological charge quantization still shows up in the limit D —cx) via the pres¬ 
ence of certain cusps in the effective potential, which are not present in the large N^. ECL 
of Ref.@]. Analogous “glued” effective potentials containing cusp singularities arise in 
supersymmetric N = 1 theories when quantization of the topological charge is imposed 
B B i- As will be discussed below, these modihcations are not essential for the local 
properties of the effective chiral potential in the vicinity of the global minimum. In this 
case, the results of Ref.0 are reproduced along with calculable l/Nc corrections. On the 
other hand, for large values of 9 and/or our results deviate from those of p. Last but 
not least, the problem of the 9 dependence in QCD is directly relevant for the construction 
of a realistic axion potential that would be compatible with the Ward identities of QCD. 
This is because an axion potential V (a) can be obtained, provided the functional form of 
the vacuum energy in QCD E^ac{9) is known, by the formal substitution 9 —> a{x)/fa- 
Our presentation is organized as follows. In Sect. 2 we describe the approach of Refs. 
0 and B to the construction of an anomalous effective Lagrangian for pure YM theory 
and QCD, respectively. Sect. 3 discusses different proposals to hnd the dimensionless 
integer-valued parameters that enter our results. A particular method suggested in 


to hx these numbers in the context of pure YM theory is presented in an appendix in a form 
adopted for the case of full QCD. In Sect. 4 we show how the heavy “glueball” degrees of 
freedom in our effective potential can be integrated out, thus yielding an effective chiral 
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potential for the light degrees of freedom. A correspondence with the VWW ECL is 
established. In Sect. 5 we discuss the vacuum properties, 9 dependence and domain wall 
solutions in the resulting effective theory. Sect. 6 is devoted to an analysis of the U(l) 
problem and properties of the pseudo-goldstone bosons at zero and non-zero 9. Sect. 7 
deals with the implications of our results for the properties of the axion and a possible 
study of the 9 dependence and new axion search experiment at RHIC. We also discuss 
the possibility of baryogenesis at the QCD scale, which seems suggestive in view of our 
results. 


2 Anomalous effective Lagrangian for QCD 


We start with recalling the construction ^ of the anomalous effective potential for pure 
YM theory (gluodynamics). It is dehned as the Legendre transform of the generating 
functional for zero momentum correlation functions of the marginal operators 
and which are hxed by the conformal anomaly in terms of the gluon condensate 

n, 0- The effective potential is a function of effective zero momentum helds h, h which 


describe the VEV’s of the composite complex helds iL, H : 


/ 


dx h 




(JdxH), 


( 1 ) 


where 


H 


2 \ ias ^yMdvr J 


H 


1 //?(«,) 2 


Wm dvr j 


( 2 ) 


and (3{as) = + > ^ym = (ll/3)Ac is the Gell-Mann - Low /3-function 

for YM theory^ , and ^ym is a generally unknown parameter which parametrizes the 
correlation function of the topological density 


lim i 

q^O . 




{0|r| 


^^GG(.)^GG(0) 


|0)yM = ivM 

da. 


YM 


( 3 ) 


The two-point function @ and other zero momentum correlation functions of a^G^, UgGG 
are dehned via the Wick type T-product by the nonperturbative part of the partition 
function \og{Z{9)/Zpy) ( Zpp stands for a perturbatively dehned partition function which 
does not depend on 6^), where 


Z{9) = Zpyex^i {—iVEy{9)} = Zpp exp 


-W(0| 


16 q!s 



( 4 ) 


by diherentiation with respect to the bare coupling constant I/Qq and 9. In Eq.(^ we 
used the fact that the vacuum energy is dehned relatively to its value in perturbation 
theory by a nonperturbative part of the conformal anomaly, for any 9. When the vacuum 

^ In what follows we will work with the one-loop /3-function. However, most of the discussion below 
can also be formulated with formally keeping the full /3-function. 
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expectation value (VEV) in Eq.(^) is defined in this way, its dependence on l/gfg is fixed 
by the dimensional transmutation formula 


(- -G"^) = const 

' Svr ' 


Mr exp 



(5) 


Here Mr is the ultraviolet cut-off mass, and the one-loop /3-function is used. It is im¬ 
portant to stress that different regularization schemes generally lead to different values of 
the constant in Eq.@ but, once specified, the VEV @ determines all zero momentum 
correlation functions of /3(as)/(4as)G^, with perturbative tails subtracted. At 6* = 0, 
Eq.(H) follows from Eq.(§) using the general relation |Q 




-G%fie) , f{e)^i-2ee^ + oie^) 


( 6 ) 


The strong assumption made in was that Eq.@ is actually covariant in 6, i.e. remains 
valid for any (at least, small) value of 6. The assumption of covariance in 6 is reproduced 
a posteriori from the effective Lagrangian, and is thus self-consistent. We note that 
covariance of Eq. (^) in 6 follows automatically within the approach suggested in . 
In fact, it is this conjectured covariance of Eq.(^ in 6 that underlies the holomorphic 
structure of the resulting effective potential (see Eq.(|^) below). Thus, we are not able 
at the moment to prove holomorphy, but instead argue that it is there based on (i) 
self-consistency of this proposal (see Sect.3), and (ii) the possibility of comparing our 
hnal formulas with the known results (such as the large W effective chiral Lagrangian 
and anomalous Ward identities in QCD, see Sect.4), and the experience from the known 
models. 

The advantage of using the combinations (@) is in the holomorphic structure of zero 
momentum correlation functions of operators G^ , GG written in terms of the H, H fields 


10 : 


lim i J dxe^'^^{0\T{H{x) H{0)}\0) = -A{H) , 

limi f dxe^'^^{0\T{H{x) H{0)}\0) = -4(iL) , (7) 

lim i f dxe^'^G0\T{H(x) 77(0)}|0) = 0 . 
g^o J 

In what follows Eqs.(|^) will be sometimes referred to as the anomalous Ward identities 
(Wrs)0. It can be seen that the n-point zero momentum correlation function of the 
operator H equals (—4)"“^(77). Multi-point correlation functions of the operator H are 
analogously expressed in terms of its vacuum expectation value (77). At the same time, 
it is easy to check that the decoupling of the fields 77 and 77 holds for arbitrary n-point 
functions of 77 , 77. This is the origin of holomorphy of an effective Lagrangian for YM 
theory, which codes information on all anomalous WI’s. 

^It should be noted that the relations (^ are not the genuine anomalous WI’s, as Eq.(^ used in (^ 
is not a Ward identity in the proper sense, see the discussion above. 
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One should note that the right hand side of the last equation in (^) does contain per¬ 
turbative contributions proportional to regular powers of However, they are irrelevant 
for our purposes, as we are only interested in the decoupling of the helds H and H at 
the level of nonperturbative effects. Holomorphy of an effective potential for 

YM theory has the same status. Thus, in contrast to the supersymmetric case where 
holomorphy is an exact property of the effective superpotential, in the present case it 
only refers to a “nonperturbative” effective potential which does not include perturbative 
effects to any hnite order in We assume that perturbative and nonperturbative effects 
can be separated, at least in principle or/and by some suitable convention HI^, As a 
result, a perturbatively dehned partition function ZpT bearing a non-holomorphic depen¬ 
dence on Qq decouples in zero momentum correlation functions. On the other hand, the 
9 dependence appears only in the nonperturbative part of Zq. Thus, the nonperturbative 
vacuum energy depends only on a single complex combination r = 1/(70 + i9bq/{327i‘^p) 
113 0- Indeed, arguments based on renormalizability, analogous to those used in Ref. 
El (see also Sect. 3), imply the relation 

(- -G‘^)e = const Re M/ exp [——5 - i-9] = const Re M/ exp [- - —r | . ( 8 ) 

Stt \ m P ) \ b ) 

(This expression coincides with the result obtained in ^ directly from the effective La- 
grangian (P).) This is exactly the origin of the relations (j^) which are obtained by dif¬ 
ferentiation of \og{Zg/Z pt) with respect to the holomorphic sources r, r. Thus, once an 
assumption of separation of perturbative and nonperturbative contributions in Eq. (P is 
made, a new complex structure emerges due to the nonperturbative origin of the 9 param¬ 
eter which combines with another parameter 1 /into the unique complex combination 
r. 

The hnal answer for the improved effective potential W{h, h) (here ’improved’ refers 
to the necessity of summation over the integers n, k in Eq.(^, see below) reads ^ 


^-iVW(h,h) 


+ 


+00 <7-1 r -y / 

E Eexp Ulo, 


n=-oo k=0 


iirV k + 


4 V 


a 


h - h 
+ ti Log 


YM 


C- 


YM , 


q 9 + 27m \ h — h 


p 


2ti 


2i 


( 9 ) 


where the constants Cymi Cym can be taken to be real and expressed in terms of the vac¬ 
uum energy in YM theory at 6 * = 0, Cym = Cym = — 2eEA^)(o) = — 2e(—6yM«s/(327r)G'^), 
and V is the 4-volume. The integer numbers p and q are relatively prime and related to 
the parameter ^ introduced in Eq.(|[) by q/p = 2^. Thus, we expect that the parameter ^ 
dehned in Eqs.(^,(|^) is a rational number. This expectation is motivated by the fact that 
it turns out to be the case in all existing proposals to hx the value of to be discussed 
in the next section, and by experience with supersymmetric models. (In all likelihood, 
irrational values of ^ would produce a non-differentiable 9 dependence for YM theory.) On 
general grounds, it follows that p = 0{Nc ), q = 0{N^). The symbol Log in Eq.(||) stands 
for the principal branch of the logarithm. The effective potential (^) produces an inhnite 
series of anomalous WI’s. By construction, it is a periodic function of the vacuum angle 
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9 The effective potential (|^) is suitable for a study of the YM vacuum as described 
above. 

The double sum over the integers n, k in Eq. (|^) appears as a resolution of an ambiguity 
of the effective potential as defined from the anomalous WTs. As was discussed in Q, 
this ambiguity is due to the fact that any particular branch of the multi-valued function 
h\og{h /corresponding to some fixed values of n,k, satisfies the anomalous WTs. 
However, without the summation over the integers n,k in Eq.(P), the effective potential 
would be multi-valued and unbounded from below. An analogous problem arises with 
the original VY effective Lagrangian. It was cured by Kovner and Shifman in ^ by a 
similar prescription of summation over all branches of the multi-valued VY superpotential. 
Moreover, the whole structure of Eq.(|^) is rather similar to that of the (amended) VY 
effective potential. Namely, it contains both the “dynamical” and “topological” parts 
(the first and the second terms in the exponent, respectively). The “dynamical” part of 
the effective potential (j^) is similar to the VY |]^ potential ~ S'log(S'/A)^'= (here S is 
an anomaly superfield), while the “topological” part is akin to the improvement 0 of 
the VY effective potential. Similarly to the supersymmetric case, the infinite sum over 
n reflects the summation over all integer topological charges in the original YM theory. 
The difference of our case from that of supersymmetric YM theory is that an effective 
potential of the form {1 /N)(j)\og[cj)/A)^ , as in the SUSY case, implies a simpler form of 
the “topological” term ~ 2mn/N {(f) — 0) with only one “topological number” n which 
specifies the particular branch of the multi-valued logarithm. In our case, we allow for 
a more general situation when the parameter Y is a rational number N = p/q. In this 
case we have two integer valued “topological numbers” n and k, specifying the branches 
of the logarithm and rational function, respectively. Our choice is related to the fact that 
some proposals to fix the values of p, q suggest that g 7^ 1, see Sect. 3. (As follows from 
Eq.(||), the values of p,q are fixed if the 9 dependence is known.) One may expect that 
the integers p and q are related to a discrete symmetry surviving the anomaly, which may 
not be directly visible in the original fundamental Lagrangian. 

It should be stressed that the improved effective potential (^ contains more infor¬ 
mation in comparison to that present in the anomalous Ward identities just due to the 
presence of the “topological” part in Eq.(^. Without this term Eq.(H) would merely be 
a kinematical reformulation of the content of anomalous Ward identities for YM theory. 
The reason is that the latter refer, as usual, to the infinite volume (thermodynamic) limit 
of the theory, where only one state of a lowest energy (for 9 fixed) survives. This state 
corresponds to one particular branch of the multi-valued effective potential in Eq.(^. 
At the same time, the very fact of multi-valuedness of the effective potential implies that 
there are other vacua which should all be taken into consideration when 9 is varied. When 
summing over the integers n, k, we keep track of all (including excited) vacua of the the- 

"‘As was explained in 6* should appear in the effective Lagrangian in the combination 9 -|- 27rn, as 
this combination arises in the original YM partition function when the topological charge quantization 
is explicitly imposed. This prescription automatically ensures periodicity in 9 with period 27r. However, 
with an additional integer k in Eq.(|^ and the above way of introducing 9 into the effective potential, 
one more invariance under the shifts 9 9 — 27(1/q arises. As will be discussed in Sect.5 in the context 

of QCD, the choice g 0 1 is not in contradiction with the known results concerning the 9 dependence for 
small values of 9. 
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ory, and simultaneously solve the problems of multi-valuedness and unboundedness from 
below of the “one-branch theory”. The most attractive feature of the proposed struc¬ 
ture of the effective potential is that the same summation over n, k reproduces the 
topological charge quantization and 2tt periodicity in 9 of the original YM theory. 

We now proceed to the generalization of Eq.(^ to the case of full QCD with Nf 
light flavors and W colors. In the effective Lagrangian approach, the light matter helds 
are described by the unitary matrix Uij corresponding to the 75 phases of the chiral 
condensate: Uij with 


U = exp 




UU^ = 1 , 


( 10 ) 


where A“ are the Gell-Mann matrices of SU{Nf), 7 r“ is the pseudoscalar octet, and = 
133 MeV. As is well known |^, the effective potential for the U field (apart from the mass 
term) is uniquely determined by the chiral anomaly, and amounts to the substitution 


e^e-iTr\ogU 


( 11 ) 


in the topological density term in the QCD Lagrangian. The rule (pTl) is valid for any 
Nf.. Note that for spatially independent vacuum fields U Eq.(pd^ results in the shift of 
0 by a constant. This fact will be used below. Furthermore, in the sense of anomalous 
conformal Ward identities []^ QCD reduces to pure YM theory when the quarks are 
“turned off” with the simultaneous substitution {G‘^)qcd {G‘^)ym and b = bqcD 
byM. Analogously, an effective Lagrangian for QCD should transform to that of pure YM 
theory when the chiral fields U are “frozen”. Its form is thus suggested by the above 
arguments and Eqs. (i),( 0 ): 

q-l 


^-iVW(h,U) _ 


n=—QO k=0 \ \ 


h Loq -h h Loq 

^ 2eE ^ 


h 


2eE. 


mV 


^ q 9 — i\og Det U + 27rn 


ft 


-hi 

^ + -VTt{MU + h.c.) 


( 12 ) 


where M = diag{mi\{^^'^^)\) and the complex fields h, h are defined as in Eq.(|^) with the 
substitution byM ^ b = (ll/3)Yc — {2/3)Nf. The integers p,q and parameter ^ {q/p = 
2^) in (pYl) are in general different from those standing in Eq.(^. Possible approaches 
to hx the values of parameters p and q in gluodynamics and QCD will be discussed 
in the next section. The constant E can be related to the gluon condensate in QCD: 
E = {bas/{S27i)G‘^), as will be clear below. We note that the “dynamical” part of the 
anomalous effective potential (|^ can be written as Wd + where 

' h DetU 


WfiiK U) = ]^hLog 


o—iO 


-Tr MU 
2 


(13) 


which is quite similar to the effective potential [0 for SQCD0. 

® For an early attempt to search for holomorphy in the effective Lagrangian framework for QCD, see 
[p^ . The problem with the approach of Ref. was that the resulting effective potential was multi-valued 
and unbounded from below. The prescription of summation over all branches of the multi-valued action 
in Eqs. (p|), (12) cures both problems. 
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Let us now check that the anomalous WI’s in QCD are reproduced from Eq.(|^). The 
anomalous chiral WI’s are automatically satisfied with the substitution (^) for any Nc, 
in accord with [Q. Further, it can be seen that the anomalous conformal WI’s of [|^ for 
zero momentum correlation functions of operator in the chiral limit niq 0 are also 
satisfied with the above choice of constant E. This is obvious from Eq.(|^, see below. 
As another important example, we calculate the topological susceptibility in QCD near 
the chiral limit from Eq.(^). For simplicity, we consider the limit of SU{Nf) isospin 
symmetry with Nf light quarks, <C Eqqd- For the vacuum energy for small 9 < ir/q 
we obtain (see Eq.(^^ below) 


E^ac{9) =-E + f cos \ + 0{m^^) 



(14) 


Differentiating this expression twice with respect to 9, we reproduce the result of |T4 


lim i f dxE'^GO\T\—GG(x) —GG(0) \ |0) = 
J f 87r Svr 


d‘^E, 


^ = T,„(f*) + oK). 


(15) 

Other known anomalous WI’s of QCD can be reproduced from Eq. ([I^ in a similar way. 
Therefore, we see that Eq.(|^ reproduces the anomalous conformal and chiral Ward 
identities of QCD and gives the correct 9 dependence for small values of 9, and in this 
sense passes the test for it to be the effective anomalous potential for QCD. Further 
arguments in favor of correctness of Eq.(^) will be given in Sect. 4, where we show that 

in the vicinity of the global minimum in 
■‘glueball” degrees of freedom, and in 


Eq.(|^ correctly reproduces the VVW ECL 
the large W limit after integrating out the heavy 
addition yields an infinite series of 1/W corrections. On the other hand, we will explain 
why we obtain a different behavior of the effective chiral potential for large values of the 
chiral condensate phases 0*. 


3 What are the values of parameters p and q ? 


In the previous section we have considered the anomalous effective potentials for YM 
theory and QCD, which involve some integer numbers p and q, with 2^ = g/p, which were 
not specified so far. The purpose of this section is to describe different proposals to fix 
the numbers p, q, which exist in the literature, and to make some comments on them. A 
related discussion can be found in the Appendix. 

Historically, the first suggestion to fix the proper holomorphic combinations of the 
fields G^ and GG was formulated in Ref. [M . For the case of pure YM theory, the authors 


proposed that fields of definite dualities dominate the vacuum, and therefore the correct 
holomorphic combinations are G^ ± iGG. Leaving aside the issue of justification of this 
hypothesis, it is of interest to discuss what values of the parameters p and q are implied 
in this scenario. As was argued in [^, the 9 dependence of the vacuum energy is fixed 


in this case by the renormalization group arguments, since for the VEV’s of interest the 
net effect of the 9 term reduces to the redehnition of the coupling constant 

1 i9 


1 




















which yields for the vacuum energy for small values of 6 

ba 


{ - 0^)0 = const Re Mq exp 

Stt 


327r2 .40' 


(17) 


which corresponds to g/p = 4/6, see Eq. (H). Thus, we see that the self-duality hypothesis 
of Ref.[|r^ implies the values p = lliV/, g = 12 , for generic odd values of Nc 7 ^ 3fc with 
some integer k. As the value of p determines the number of different non-degenerate 
vacua in the theory [^, we end up with llA"c vacua, which may look strange. This is in 
contrast to the case of supersymmetric YM theory where the holomorphic combinations 
are known to be ± iCG, but the number of vacua is for any number of colors. 
This may be understood in terms of the renormalization group arguments similar to Eqs. 
(|^),(^) (with the substitution (G^) —(AA)) as a result of the interplay between the 
integer valued /3-function bsvM = 3Yc, which is determined by the zero modes alone and 
has a geometrical meaning, and the dimension d = 3 of the gluino condensate. It appears 
that this conspiracy is very specific to supersymmetric theories. It is interesting to note 
in this reference that if for some reason only the zero mode contribution b = 4Yc, instead 
of the full 6 = (4 — l/3)Ac, were to be retained in the /3-function, Eq.([l7D would imply 
Nc vacua. However, we are unable at the moment to see any compelling reason why 
such substitution should be made. Thus, it remains unclear whether or not the appealing 
choice G^±iGG , p = Nc can be compatible with the renormalization group and conformal 
anomaly for non-supersymmetric YM theory. 

Another approach to the problem of the number of vacua and proper holomorphic 
combinations of the fields G^ and GG is based on the analysis of softly broken SUSY 
theories |T^ , which is under theoretical control as long as the gluino mass is much smaller 


than the dynamical mass scale: -C Asym- A rather detailed discussion of this scenario 


has been recently given by Shifman using supersymmetric gluodynamics with = 3 
as an example. In the limit of small rUg the VEV of the holomorphic combination G^+iGG 
is proportional to the VEV mg{XX) where the gluino condensate (AA) is to be calculated 
in the supersymmetric limit rUg = 0. The 6 dependence of the latter is known 0: 
(AA) ~ exp{i9/Nc + 27ik/Nc) ,k = 0,1,Nc — I, which corresponds to Nc degenerate 
vacua. When nig ^ 0, the vacuum degeneracy is lifted. For Nc = 3 and 0 = 0, we have 
one state with negative energy E = —mgAgyM^ degenerate states with positive 

energy E = {l/2)mgA^YM- The former is the true vacuum state of softly broken SUSY 
gluodynamics, while the latter are metastable states with broken GP. The lifetime of the 
metastable states is large for small nig, and decreases as nig approaches Asym ■ When 0 
is varied, the three states intertwine, thus restoring the physical 27r periodicity in 0. This 
picture suggests the values p = Nc, g = 1. 

The problem with the above SUSY-motivated scenario is that the genuine case of pure 
YM theory corresponds to the limit nig 3> Asym which is not controlled in this approach. 
Moreover, the conformal anomaly in softly broken SUSY gluodynamics is different from 
that of pure YM theory. On the other hand, it is clear from the above discussion that the 
conformal anomaly and dimensional transmutation are very essential for the analysis of 
the 0 dependence in gluodynamics. Perhaps, it is worthwhile to mention that the value 
p = llVc follows also within a non-standard non-soft SUSY breaking suggested recently 
^ as a toy model to match the conformal anomaly of non-supersymmetric YM theory 
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at the effective Lagrangian level. 

With these reservations, it is nevertheless reasonable to expect that the above SUSY- 
motivated scenario is close to what actnally happens in the deconpling limit nig ^ Asym- 
Two different versions of this scenario may be expected. First, it may happen that at all 
generic valnes of 9 there exist one trne vacnnm of lowest energy pins (p — 1) metastable 
CP-violating vacna, which are separated by potential barriers and intertwine when 6 
evolves. Another possibility is that metastable vacna exist only in the vicinity of a level 
crossing point in 6, while for other valnes of 6 they become the saddle points or maxima 

|T^. In one of these forms, snch a pictnre seems to be needed to match the Witten- 


Veneziano resolntion of the U(l) problem. The scenario discnssed in implies that 
the nnmber of vacna p remains W in the limit nig S> Asym, bnt it is conceivable that 
an additional level splitting occnrs with passing the region nig ~ Asym where the SUSY 
methods become inapplicable. Actnally, the pictnre arising in onr approach will be jnst 
in this vein (see Sect.5). As will be discnssed there, which of the above two versions is 
realized is mostly determined by the valne of parameter q. When q ^ 1, metastable vacna 
exist for all valnes of 6, as it happens in the SUSY scenario. Fnrthermore, it will be 
argned that large lifetimes of metastable states, necessary for this scenario to work, are 


ensnred parametrically - a fact which is not seen []Ig in the SUSY-motivated pictnre. 

One more possible approach we wish to discnss is based on an idea formnlated some 
time ago by Kiihn and Zakharov (KZ) These anthors have snggested that in QCD 
with massless qnarks nonpertnrbative matrix elements shonld be holomorphic in the Panli- 
Villars fermion mass Mp^. Assnming this kind of holomorphy, they have proposed to relate 
the proton matrix element of the topological density {p\GG\p) to the matrix element 
{p\G‘^\p) which is fixed by the conformal anomaly. However, it is not easy to separate 
pertnrbative and nonpertnrbative contribntions to the latter dne to a non-trivial proton 
wave fnnction. This may be a potentially problematic point when the KZ holomorphy is 
considered for the matrix elements]^. On the other hand, for vacnnm condensates and zero 
momentnm correlation fnnctions the nonpertnrbative contribntions can be systematically 
singled ont, at least formally |T^. In this respect, the latter objects are simpler than 


the hadron matrix elements, and thns appear preferable for testing the KZ holomorphy. 
This issne was addressed in where a method similar to that of Ref. [pl[] was nsed 


in the context of pnre YM theory to relate the zero momentnm two-point fnnction of 
GG to that of G^. In this approach, pnre glnodynamics was considered as a low energy 
limit of a theory inclnding a heavy qnark, while holomorphy in the physical fermion 
mass m —> cx) was argned to hold basing on deconpling argnments. Appealing to this 
holomorphy, it was snggested that for the case of pnre YM theory the parameters of 
interest are p = SbYM = and q = 8, for odd Nc- Fnrthermore, as the anomalous 

conformal WFs for the operator G^ |j^ are covariant in 6, one can conclude that the 
Kiihn-Zakharov holomorphy, if it holds, indeed implies covariance of Eq.(|^) in 6, and thus 
leads to the holomorphic effective potential for glnodynamics, Eq.(^. 

An inverse route was undertaken in [^. In this paper, the starting point was the 
holomorphic effective potential (^) for pure YM theory, with unspecified parameters p 
and q. Again, the idea was that the meaning of this holomorphy can be clarified by 


“V.I. Zakharov, private communication. 
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coupling pure YM theory to a very heavy fermion with mass m ^ oo, while the values of 
parameters p and q would be hxed in this case by some kind of consistency conditions. In 
contrast to the previous approach, it was suggested in m to introduce a heavy fermion 
directly at the effective Lagrangian level by using the “integrating in” procedure familiar 
in the context of SUSY theories [^, ^]. Thus, the “integrating in” method is used here to 
construct an effective Lagrangian for the system (YM + heavy fermion) from the effective 
Lagrangian for pure YM theory. The consistency condition suggested in is that the 


holomorphic structure of the latter should arise from the holomorphic structure of the 
former, assuming the standard form of the fermion mass term. Then, the “integrating 
in” method plus the above consistency condition are found to select the only possible 
values p = SbyM = H-Yc and q = 8. These results coincide with the ones obtained with 
the approach of [^. Thus, the “integrating in” method seems to suggest the effective 
Lagrangian realization of the Kiihn-Zakharov holomorphy. The agreement of these two 
lines of reasoning is encouraging, and shows that different assumptions made in [|^] and 
are at least consistent with each other. 

Finally, we would like to comment on another related development. Very recently, Wit¬ 


ten 1^ has shown how the qualitative features of the 9 dependence in non-supersymmetric 
YM theory - such as a multiplicity of vacua ~ W, existence of domain walls and exact 
vacuum doubling at some special values of 6^ - can be understood using the AdS/CFT 
duality. The latter provides a continuum version of the strong coupling limit, with 


a hxed ultraviolet cutoff, for YM theory with W — »• cxo, —>■ cxd. As was shown in 


241, in this regime the 9 dependence of the vacuum energy in YM theory takes the form 


E^ac{,9) = C mm{9 + 2nkY + 0(1/^,) 

k 


(18) 


where C is some constant. We would like to make two comments on a comparison of our 
results with the picture advocated by Witten in the large Yc limit. First, we note that 
the structure of Eq.([T8|) agrees with our modihed dehnition of the path integral including 
summation over all branches of a multi-valued (effective) action. Indeed, Eq. (P!B| ) suggests 
the correspondence 


C min {9 -|- 27rfc)^ lim 

k V^oo 




^ g-VC(6»+27rfe)2 

L k 


(19) 


using the dehnition of the vacuum energy through the thermodynamic limit of the path 
integral. With this dehnition which prescribes the way the volume V appears in the 
formula for the vacuum energy, the correspondence (0) appears to be the only possible 
one. On the other hand, the latter expression has exactly the structure that arises with our 
dehnition of the improved ehective potential (|^). Therefore, our prescription of summation 
over all branches of a multi-valued (ehective) action seems to be consistent with the picture 
developed by Witten using an approach based on the AdS/CFT correspondence. In 
particular, our picture of bubbles of metastable vacua bounded by domain walls considered 
in the context of QCD in Sect. 5 is in qualitative agreement with that suggested by Witten 
^1 for the pure YM case. 

Second, one may wonder whether the approach of Ref. |2^ can provide an alternative 
way to hx the parameters p, q of interest. We note that Eq.([T8|) indicates a non-analyticity 
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at the values 6c = 7i {mod27i) only, where CP is broken spontaneously. If the technique 
based on the AdS/CFT duality could be smoothly continued to the weak coupling regime 
of non-supersymmetric YM theory, this would result in the values g = l,p ~ Nc- However, 
the possibility of such extrapolation is unclear, as for small A = g^M^c the background 
geometry develops a singular behavior and the supergravity approach breaks down. There 
might well be a phase transition ||^ when the effective YM coupling Qym is reduced. 


That such a phase transition should occur in the supergravity approach to QCD^ was 
argued in . Other reservations about the use of the supergravity approach to the non- 
supersymmetric YM theory in D=4 have been expressed in where no perturbative 


indication was found for decoupling of unwanted massive Kaluza-Klein states of string 
theory. On the other hand, there exists some evidence from lattice simulations that 
a critical value of 6 moves from = vr in the strong coupling regime to 6c < n in 
the weak coupling regime |29]. In terms of parameters p, g, such a case corresponds to 


g 7 ^ 1. Therefore, we conclude that if no phase transition existed in the supergravity 
approach, the results of Refs. [1^, would be in conflict with the latter which would 
imply p = 0{Nc),q = 1. In this case, the assumptions made in |T2|, would have to 


be reconsidered. Alternatively, there might be no conflict between the two approaches if 
such a phase transition does occur. 

To summarize, at the qualitative level we expect a vacuum structure similar to that 
suggested by the SUSY scenario. As for the quantitative results for the parameters p 
and g, different lines of reasoning lead to generally different answers. The self-duality 
hypothesis, Kiihn-Zakharov - type arguments and non-standard SUSY breaking toy model 
all suggest p = SbyM = llA"c with g = 8 or 12 for the pure YM case. The compatibility 
of the appealing choice ± iGG , p = Nc , q = I, suggested by the soft SUSY breaking 
scenario, with the renormalization group and conformal anomaly for non-supersymmetric 
YM theory remains unclear. It is conceivable that our current understanding of the 
effective Lagrangian is incomplete, and a more careful analysis - perhaps along the lines 
of Refs. [^, Q - will solve the puzzle (if it is a puzzle) of “extra 1/3 ”, thus favoring the 
SUSY-type scenario. Alternatively, the numbers p = llYc , g = 8 (for odd Nc) suggested 
by the methods of may be the correct answer, though perhaps somewhat “counter¬ 

intuitive” . For these reasons, below we will keep the general notation p, g, while a separate 
analysis will be given in cases when the concrete values of p, q are essential. The reader 
is referred to the Appendix for the details of the “integrating in” method of [|^ adopted 
to the case of full QCD, which suggests the values p = 3b = llYc — 2Nf , g = 8 . 


4 Effective chiral Lagrangian for finite Nc 

The anomalous effective potential ( [T^ ) contains both the light chiral fields U and heavy 
“glueball” fields h , h, and is thus not an effective potential in Wilsonian sense. On the 
other hand, only the light degrees of freedom, described by the fields U, are relevant for 
the low energy physics. An effective potential for the U , U fields can be obtained by 
integrating out the h, h fields in Eq.(|T^). It corresponds to a potential part of a low 
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energy Wilsonian effective Lagrangian for energies less than the glueball massesf]. The 


transition from the effective potential ([T^) for the U , h helds to a Wilsonian effective 
potential for the U fields by integrating out the h , h fields is analogous to the transition 
^ from TVY effective Lagrangian 0 for SUSY QCD to the Affleck-Dine-Seiberg 


low energy effective Lagrangian. The purpose of this section is to obtain such an effective 
potential for the light U, U helds by integrating out the h , h helds in Eq.(|^. 

To hnd this ehective potential for the light U , U helds, we make two observations. 
First, the mass term in Eq.(|^ does not couple to the “glueball” helds, and is thus 
unessential for integrating them out. To simplify the subsequent formulas, in this section 
we will omit the mass term, and add it at the end of calculation. Second, the only 
remaining ehect of the light matter helds U, U is formally reduced, as was discussed in 
Sect.2, to the redehnition (shift) ( PT]) of the 6 parameter, and the changes of the numerical 
parameters in the ehective potential (^) for pure YM theory. Therefore, for space -time 
independent helds U , U we can integrate out the “glueball” helds h, h in Eq. ( [T^) in the 
same way as the 6 dependence of the vacuum energy for the pure YM case was found in 0 . 
For the sake of completeness, this calculation will be repeated below for the present case 
of QCD. As before, we will keep the total space-time 4-volume hnite, while a transition 
to the thermodynamic limit U —> cxd will be performed at the very end. 

We start with introducing the “physical” real helds p, uj dehned by the relations 


h = 2EeP+^^ , h = 2Et 




( 20 ) 


(This dehnition implies lU(a; -|- 27r) = W{uj). As will be seen, this condition of single¬ 
valuedness of the ui held is satished with the substitution (p^).) Then, for the “dynamical” 
part of Eq.(]^) we obtain 


— ^^(hLog—— + h Log —^ j = —iUE [(p — 1) coso; — a; sino;] . (21) 

The summation over the integers n in Eq. (p!2[) enforces the quantization rule due to the 
Poisson formula 




sm 00 — n 


( 22 ) 


which rehects quantization of the topological charge in the original theory. Therefore, 
when the constraint (|2^) is imposed, Eq.(pT|) can be written as 


Using 


xV / h — h \ p 

— h Loq —— -|- h Loq —— = —iVEe^ (p — 1) coso; + in-oo . 
4 \ ^2eE ^2eE yy ) ^ 


we 


put Eq. ([I^) (with the mass term omitted) in the form 


+00 g-l 

^-ivw ^ ^ ^ (5(UE - e^sino;—n) exp 

n=-oo k=0 P 


P 


-iVE e^{p — 1) cos 00 + in \ 6k + - 00 


(23) 


(24) 


^ Such an object is not precisely Wilsonian effective action in the usual sense as it does not involve 
e.g. the vector p,oj mesons whose masses are compatible with that of the g'. 
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where we denoted 


6k = 6 — i log Det U + 27r- k . 

q 


(25) 


To resolve the constraint imposed by the presence of 5-fnnction in Eqs. (P^ , (|^) , we 
introdnce the new held $ by the formula 


6{VE - sin a; — n) oc D ^ exp [ i$l/Ee^sina; — - n | 

p V q J 

Going over to Euclidean space by the substitution iV —> V, we obtain from Eqs. 

^ ( +00 q -1 

hE(t/, p, o;, <h) = —— log < ^ exp \—VEe^{{p — 1 ) cos a; — $sina;} 

^ [_n=—oofe=0 

+ in (9 - ilogDetU + 27ik- + - u - - <!>] -e^ 

\ q q q J ee 


(26) 


(27) 


Here we introduced the last term to regularize the inhnite sum over the integers n. The 
limit e —0 will be carried out at the end, but before taking the thermodynamic limit 
E —>• cx). Note that Eq.fpTl) satishes the condition iy(a; + 27r) = W{ui) which should hold 
as long as uj is an angle variable. We also note that the periodicity in 9 with period 2ti is 
explicit in Eq.(^. 

To discuss the thermodynamic limit E —> cx) we use the identity 


03(1/, x) = 


+ 00 

Z! exp 


k=—oo 


{v + kf 


X 


+ O0 


— y) exp[—+ 2ilvT:] 


(28) 


l = — OQ 


and transform Eq.(|27|) into its dual form 

1 ( +00 q-1 

fE(f/, p, cn, <h) = —— log < Z Z 6xp [—EEe^ |(p — 1) cos a; — <hsina;| 

^ [_n=—cxDfc=0 

- — (9 - ilogDetU + 2nk^ + ^ UJ - ^ ^ - 27m] 

Ae \ q q q J 


(29) 


where we have omitted an irrelevant inhnite factor ~ in front of the sum. Eq.(M) 
is the hnal form of the improved ehective potential W, which is suitable for integrating 
out the “glueball” helds p, uj, along with the auxiliary held <h. To this end, the function 
W should be minimized in respect to the three variables p,uj and $, with U hxed. In 
spite of the frightening form of this function, its extrema can be readily found using the 
following simple trick. As at the extremal points all partial derivatives of the function W 
vanish, we hrst consider their linear combination in which the sum over n, k cancels out. 
We thus arrive at the equations 


dW 

dp 

dW dW 


= Ee^{pcosui — ^smuj) = 0 


+ ^—= —Ee^(psina; + "hcoso;) = 0 , 


duj 9$ 


(30) 
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which is equivalent to = 0. (We do not consider here the case p ^ —oo which 

would also solve Eqs. (|30D , see [|T^.) Therefore, these equations have the only solution 


(p) = 0 , ($) = 0 , 


(31) 


while the minimum value of the angular held oj is left arbitrary by them. The latter can 
now be found from either of the constraints dW/duj = 0 or dW/d^ = 0, which become 
identical for (p) = (<h) = 0. The resulting equation reads 


+00 g —1 


E E( 9 — i log Det U + 2'Kk - — 2Trn + -ui + 2e- 


P 


sino; 


n=-oo k=0 


q 


q 


p 

2 ' 


X exp <VE cos uj 


9 — i log Det U + 2'Kk - + - a; — 27rn ) }■ = 0 , (32) 


As 


in which we have to take the limit e —>■ 0 at a hxed 4-volume V. One can see that 
non-trivial solutions of Eq. (^) at e —>■ 0 are given by 

(u}i = --(9 — Hog Det U)-i -Z-|-27rr , Z = 0,1,... ,p — 1 ; r = 0,±l,... (33) 

p p 

Eq. (P5D shows that p physically distinct solutions of the equation of motion for the uj held, 
while the series over the integers r in Eq. (|33D simply rehects the angular character of the 
OJ variable, and is thus irrelevant. Substituting Eq. (|33D back to Eq.(^9D and restoring the 
mass term for the U held, we obtain the ehective potential for the light chiral helds U , U 

B- 




V^oo V 


( p-1 

- 

S Eexp 

VE cos 

{l=0 

- 


--(9 — i log Det U) + — l 

p p 


-VTr{MU + M+U+) 


(34) 


Eq. ( p4D is the hnal form of the ehective potential for the chiral held which is valid for any 
value of 9. As was mentioned in the beginning of this section, this form of the ehective 
potential (|5^ could be read oh the formula ^ for the vacuum energy as a function of 9 
in pure YM theory, with the substitution of the parameters p, q by their values in QCD, 
adding the mass term for the chiral held, and making the shift (^) of the 9 parameter. 
The meaning of the resulting expression is diherent, however. What was the the vacuum 
energy as a function of 9 describing p diherent vacua in YM theory becomes the ehective 
potential for the light chiral helds. At hrst sight, it could be expected that the resulting 
ehective chiral Lagrangian has the same number p of diherent vacua. It turns out that this 
naive expectation is wrong: the number of diherent vacua in QCD is determined by the 
number of havors Yj, at least as long as Nf < N^. The reason why the number of vacua 
diher from the number of branches of the ehective potential is the angular character of 
the corresponding chiral degrees of freedom. We hope that the last sentence will become 
clearer in the next section when we consider the concrete examples. 
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The peculiarity of the resulting effective potential (|51D is that it is impossible to 
represent it by a single analytic function by directly performing the limit 1 / —> oo in (|3^ . 
In the thermodynamic limit Id —cx) the only surviving term in the sum in Eq.(|S^) is 
the one maximizing the cosine fnnction. Thus, the thermodynamic limit selects, for given 
value of 6 ^ — ilogDetU, a corresponding value of I, i.e. one particular branch in Eg. (|3^ . 
The branch structure of Eg . (13^ shows up in the limit E —> cxo by the presence of cusp 
singularities at certain values of 6 — i log Det U. These cusp singularities are analogous to 
the ones arising in the case of pure gluodynamics for the vacuum energy as function 
of 6, showing non-analyticity of the 6 dependence at certain values of 6. In the present 
case, the effective potential for the light chiral helds analogously becomes non-analytic at 
some values of the helds. The origin of this non-analyticity is the same as in the pure 
YM case - it appears when the topological charge quantization is imposed explicitly at 
the effective Lagrangian level. 

The general analysis of the effective potential (|5^) will be given in the next section, 
while here we consider the case when the combination 6 — i log Det U is small, and thus 
the term with I = 0 dominates. We obtain for this case 


W^ff\u,U+) = -Ecos 


--(6 — i log Det U) 
P 


--Tr {MU + M+U+) . (35) 


Expanding the cosine (this corresponds to the expansion in q/p ~ 1/iVc), we recover 
exactly the ECL of 0 at lowest order in l/W (but only for small 6 — i log Det U < Ti/q), 
together with the “cosmological” term —E = —{bas/{32n)G^) reqnired by the conformal 
anomaly: 

W^i=f^\u,U^) = -E-^^^^^{e-ilogDetUf--Tr{MU + M+U+) + ... , (36) 

where we used the fact that, according to Eq.(^, at large W E{q/p)‘^ = —{v^)ym where 
< 0 is the topological susceptibility in pure YM theory. Corrections in 1/iVc 
stemming from Eq.(|35|) constitute a new result. Thus, in the large W limit the effective 
chiral potential (^) coincides with that of @] in the vicinity of the global minimum. At 
the same time, terms with / 7 ^ 0 in Eq.(|3^) result in different global properties of the 
effective chiral potential in both cases q = f and g 7 ^ 1 in comparison with the one of 
Ref. 0, see below. 


5 Theta dependence, metastable vacua and domain 
walls 

In this section we analyse the picture of the physical 9 dependence and vacuum structure 
stemming from the effective potential (|3^ . This is where we encounter the main difference 
of our results from the scenario of ^ . The origin of this difference is the branch structure 
of the effective potential (^^, with the prescription of summation over all branches. As we 
have mentioned earlier, this effective potential has cusp singularities at certain values of 
the fields, whose origin is the topological charge quantization in the effective Lagrangian 
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framework. It is therefore clear that these cusp singularities can not be seen in the usual 
treatment of the effective chiral Lagrangian, which deals from the very beginning with 
quark degrees of freedom only without imposing quantization of the topological charge. 
Thus, the different global form of the effective chiral potential as a function of the chiral 
condensate phases 0 *, with cusp singularities at certain values of the phases 0 *, is the 
hrst essential difference of our picture from that of Ref.[^. Another important difference 
appears when q ^ 1. In this case, we will hnd metastable vacua, separated by the barriers 
from the true physical vacuum of lowest energy, in the whole range of variation of 6. 
Properties of these metastable vacua will be discussed below. Furthermore, for the same 
case g 7 ^ 1 we will argue that the vacuum doubling at the points 9 = {2k + l)vr/g occurs 
irrespective of the particular values of the light quark masses (no Dashen’s constraint, see 
below). 

It is convenient to describe the non-analytic effective chiral potential (^) by a set 
of analytic functions dehned on different intervals of the combination 9 — i log Det U. 
Thus, according to Eq. (|3^ , the inhnite volume limit of the effective potential for the 
helds U = diag{expi(f)g) is dominated by its /-th branch: 

IFjil =-Ecos (-- 6 ^ ++ —/) - ^Mjcos0j , / = 0,1,... ,p - 1 (37) 

jj \ p p ^ p j ^ 


if 


{21 


1)- < 9 

q 


< ( 2 / + 1 ) 


TT 


q 


(38) 


This can be viewed as the set of “p” different effective potentials describing different 
branches in Eq.(|T^. The periodicity in 9 is realized on the set of potentials (|37D as a 
whole, precisely as it occurs in the pure gauge case @] where different branches undergo 
a cyclic permutation under the shift 6 * —> 6 * — 27r. As is seen from Eq. (|37D , the shift 
9 ^ 9 — 271 transforms the branch with I = k into the branch with I = k + q. In addition, 
as long as q ^ 1 , there exists another series of cyclic permutations corresponding to I = k 
and I = k + 1 in the above set, which are related to each other by the shift 9 —>■ 9 — 27i/q. If 
the number q were g = 1, the two series would be, of course, the same. As was mentioned 
in Sect.2, a value g 7 ^ 1 thus implies some discrete symmetry arising at the quantum 
level. (Although the periodicity in 9 with period 27r/g, g 7 ^ 1, rather than 27r may look 
surprising, it is not the hrst time we encounter such a situation. We would like to note 
in A^ = 2 Seiberg-Witten theories with quarks the 9 dependence for Nc = 2 has period 
TT, and not the “standard” 27r |pT|. A similar behavior was argued to hold in some 2D 


models on the lattice p^.) In what follows we will discuss both cases g = 1 and g 7 ^ 1. 
Consider the equation of motion for the lowest branch / = 0 : 


sin(-9--^(j)}i=-^sm(j)i , i = l,...,Nf (39) 

\p P / g A 

with the constraint (|38D with Z = 0 . At lowest order in 1 /W this equation coincides with 
that of 0. For general values of Mi/E, it is not possible to solve Eq. (|3^ analytically. 
However, in the realistic case -C 1, ~ 1 where = {p/q)Mi/ E, the approx¬ 

imate solution can be found. Neglecting the 0(eu,£d) terms in the phases 0* (in this 
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approximation we deal with the case Nf = 2), we obtain 


^(- 0 ) 

Eu sin 


0 

e 

■ Ai=o) 

£dSin0^ 


( 40 ) 


One sees that the constraint (|38[) is automatically satisfied. The solution of Eqs. (|40|) reads 

7 = 0 ) _ _ sin 9 


sm I 


sm < 


[m^ + + 2mumd cos 9] 


,(«= 0 ) 


rUu sin 9 


[m^ + + 2mumd cos 

sin0^^=°^ = 0{eu,£d) ■ 


+ 0{eu, Ed) ; 
+ 0{Eu, Ed) , 


(41) 


Thus, the solution for the I = 0 branch coincides with the one of Ref. P] to the leading 
order in £u,£d- Let us now concentrate on the case q = 1- For the next I = 1 branch we 
obtain, instead of (^Of ), 


01 '-) = 0 

+ = 9-27, (42) 

Eu sin sin 

One can easily see that the solution of these equations can be obtained from the 

previous one: ~ Obviously, solutions for branches with I > 1 will 

coincide with one of the two solutions modulo 27r. Furthermore, while the 

first solution (pu~d^'^ defines the location of the global minimum of the effective potential, 

the additional solution (pu~d^'^ saddle point, see Fig. |l| for the form of the effective 

potential at 0 = 0 for Nf = 2 and q = 1. Such a saddle point of the effective chiral 
potential may be of importance for cosmology and/or the physics of heavy ion collisions. 
We will not discuss these issues in this paper, but hope to return to them elsewhere. 

We have checked numerically that in the case q = 1, Nf = 3 with physical values of 
the quark masses we stay with one physical vacuum at all values of 9. In particular, no 
(stable or metastable) vacua appear at 0 = tt. Thus, for g = 1 the counting of vacua in our 
approach agrees with that of Ref.p]. On the other hand, when the Dashen’s constraint 

m 

mumd > ms\md - ruu] (43) 


is satisfied, metastable vacua appear in the vicinity of the point 9 = 7i. In particular, in 
the SU{2) limit ^ rUg the metastable vacuum exists in a very narrow region 

near 9 = ii. When 9 becomes exactly equal tt, the two vacua are exactly degenerate. This 
agrees with the picture of Ref. 0 . 

Finally, we consider the case g = 1 with Nf = 3 light flavors of equal masses. In this 
case, the metastable vacua exist in the extended region of 9 from 7r/2 to 3tt/2, analogously 
to what was recently found by Smilga [0 for the VVW potential in the same limit. The 
resulting picture of the 9 dependence of the vacuum energy is shown in Fig. 0. 
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Figure 1: Two flavor effective potential for q = 1. Indicated points are a) minimum and 
b) saddle point. 



Figure 2: Vacuum energy as a function of vacuum angle 6 for the case g = 1 with Nf = 3 
light flavors of equal masses. The solid and dashed lines represent minima and the dotted 
line is the intermediate saddle point. 
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Let us now consider the case q ^ 1. For the simplest situation of isospin SU{Nf) 
symmetry with masses rrii = m ^qcd, the lowest energy state is described by 


di=o) 


= -L , = 

I\f 


-E 


MNf cos + 0(mj) > ^ ^ ; (44) 
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Nf qNf 
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27r 


Nf qNf^ 


+ 0(m' 


Tl ^ 371 

- <9 < — 


q 


q 


(45) 

etc. Thus, the solution coincides with the one obtained by VVW at small 9 <77 jq 
up to 0{m‘^) terms. However, at larger values of 9 the true vacuum switches from (|4^) 
to (]45|) with a cusp singularity developing at 9 = vr/g. Here we remind the reader that 
the “standard” location of the hrst critical point 9 = tt corresponds to the particular case 
g = 1 in our general formulas. On the contrary, in the scenario @] the solution remains. 


if g 7 ^ 1, analytic at this point, and for vr/g < 9 < Svr/g has an energy larger than in (^Sf) . 
On the other hand, if g = 1, the picture of Ref.|P] is reproduced. Moreover, the number 
of different solutions (which may or may not be metastable vacua, depending on the signs 
of second derivatives of the potential) is precisely Nf (or qNf if g 7 ^ 1) as the phases are 
dehned modulo 27r, and thus only the hrst Nf {qNf) terms in the series (^ , ^5]) become 
operative. 

The interesting feature of the case g 7 ^ 1 is that the vacuum doubling at the points 


9k — (2A: + 1) 


71 


q 


fc = 0,1,... ,p -1 


(46) 


holds irrespective of the values of the light quark masses. This can be seen from the fact 
that the equations of motion for any two branches with I = k and I = k + 1 from the set 
(p7|) are related by the shift 9^9 — 271/q. Thus, the extreme sensitivity of the theory to 
the values of the light quark masses in the vicinity of the critical point in 9 is avoided in 
our scenario if g 7 ^ 1 , while the location of the critical point is given by 9c = 71 /q instead 
of the “standard” 9c = tt . (A similar situation was argued to hold in 2D CP^ models on 
the lattice in the weak coupling limit |^.) Another interesting feature of the scenario 
g 7 ^ 1 is the appearance of metastable vacua which exist for any value of 9, including 
9 = 0. For the physical values of the quark masses, we hnd g — 1 additional local minima 
of the effective chiral potential, which are separated by barriers from the true physical 
vacuum of lowest energy. For the illustrative purpose, we present in Fig. |] the effective 
potential at 6 * = 0 for g = 2 , A^/ = 2 . 

The existence of additional local minima of an effective potential for the case g 7 ^ 1 
leads to the well known phenomenon of the false vacuum decay |^. For the effective chiral 
potential (|3^ this effect and its possible consequences in the context of axion physics were 
briefly discussed in [1^. Below we present a somewhat more detailed discussion of this 
issue. 

We discuss the problem in a simplihed setting by considering the isospin SU{Nf) 
limit with equal (and small) fermion masses and taking all chiral phases equal, (pi = 0 , 
i.e. restricting our analysis to the “radial” motion in the (/)-space. In this setting, the 
problem becomes tractable in the spirit of Ref. ||33||. In what follows, we only consider 
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Figure 3: Effective potential for q=2 and Nf = 2. The point indicated as a) is the global 
minimum. The two points indicated as b) are identihed. 


transitions between a metastable state of lowest energy and the vacuum. To calculate 
the wall surface tension a, it is convenient to shift the vacuum energy by an overall 
constant such that the metastable state has zero energy, and to rescale and shift the 
chiral held (j) —>■ {2/— 7r/{qNf) in order to have the standard normalization of 
the kinetic term and symmetrized form of the potential. With these conventions, the 
effective potential for 6 = 0 becomes 





1 — cos 
1 — cos 





/W = Nf 



- Mf{(()) if 0 > 0 

- Mf{(j)) if 0 < 0 

TT \ / 27r 

- — cos - 

qNfj \qNf 


(47) 


The effective potential (|47|) has a global minimum at 0 = nfff/{2q^JWf) and a local 

minimum at 0 = —7rfTf/{2q^Jl^), with a cusp singularity between them (see Fig. 

We note that analogous “glued” potentials were discussed for SUSY models in a similar 
context 1^ ^ 3- 

A few comments on the above effective potential are in order. First, we note that the 
potential barrier is high (~ (G^) ) and wide, while the energy splitting for 6 = 0 


AE 


rriqNf 


(TT) 






(48) 


21 





















<l> 


Figure 4: Effective potential at 0 = 0 for equal chiral phases (j) = (l)u = 4>d = 4>s- 


is numerically small in comparison to the leading term proportional to the gluon con¬ 
densate. Therefore, the thin wall approximation of Ref.||^ is justihed for the physical 
values of the parameters that enter the effective potential. Indeed, it is easy to check 
that the necessary condition Sa/i ^ AE (where a is the wall surface tension given 
by Eq.(|^ below, and fi is the width of the wall) is satisfied to a good accuracy for both 
alternative choices p = llA'), — 2Nf = 27, g = 8 or p = iVc = 3 , g = 1. (In the latter case, 
a metastable state does not exist for Nf = 3 and 9 = 0, but is possible for, say, Nf = 5.) 
Second, we would like to comment on the meaning of the cusp singularity of the effective 
potential (^). As was mentioned earlier, the cusp arises as a result of integrating out the 
“glueball” degrees of freedom, which were carrying information on the topological charge 
quantization. For an analogous situation in the supersymmetric case, it was argued 0 
that the cusp, where the adiabatic approximation breaks down, provides a leading con¬ 
tribution to the wall surface tension. In our case, we expect the difference of the surface 
tensions for the potential (|47| ) and a potential where the cusp is smoothed to be down 
by powers of N^. The reason is that the domain wall to be discussed below is in fact the 
p' wall, while on the other hand, a coupling of the p' to the glueball fields near the cusp 
would yield the above suppression. 

Explicitly, the domain wall solution corresponding to the effective potential ( 
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, X > Xo , (49) 
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where xq is the position of the center of the domain wall and 





(50) 


is the width of the wall, which turns out to be exactly equal to the rj' mass in the chiral 
limit, see Eq. ( 0 ) below. This suggests the interpretation of the domain wall P5|) as the 
Tj' domain wall. The solution as a function oi x — xq is shown in Fig. (^). Its hrst 
derivative is continuous at a: = a:o, but the second derivative exhibits a jump. 

The wall surface tension can be easily calculated from Eq. (|T7|) . For 6* = 0 it is 


a = 




q^Nf V 327r 
Eq.(pD should be compared with the formula 


cos^j + 0 (mj 2 ). 


(51) 


TT 


a = 


372(1- ^ 

V 373, 




(52) 


found by Smilga |]^ for the the wall surface tension at 6 ^ ~ tt for the VVW potential 
with Nf = 3 with equal fermion masses. A distinct difference between these two cases 
is the absence of the chiral suppression ~ rUg in Eq. (|5TD , which apparently would make 
penetration through the barrier even more difficult in comparison to the VVW potential. 
Another difference is the large Nc behavior of Eqs. (^) and (^). For hxed rUg ^ 0, the 


surface tension scales as 0 {Nc) for both cases, while in the limit rUg —> 0 hrst and then 


Nc —>■ oo the surface tension vanishes for Eq. 


and scales as for Eq. (^). 
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The quasiclassical formula for the decay rate per unit time per unit volume is 


^ , 277rV\_ , 

1 oc exp ( ) = exp(-h 4 ). 


( 53 ) 


^ 2(A£)3^ 

Using (1^) , one hnds that in the VVW scenario the lifetime of the metastable state at 
zero temperature is much larger than the age of the Universe 
very different result^ 


19]. In our case, we find a 


^4 = 


3V 2U vrV 


(l-cos 

q^Nj 

M3 



27 TrViV; 
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m„ 


(TT) 


(54) 


Eq. 


shows that the false vacuum decay is suppressed parametrically by a factor ~ 
i^QCD/fnq)^, which should be compared to a factor ~ {-^sym/^qY in the soft SUSY 
breaking scheme While the latter ceases to yield a suppression with approaching 

the decoupling limit mg ^ Asym, the former is a real suppression factor for QCD. As 
should be expected, it tends to infinity (i.e. vacua become stable) when goes to zero. 
On the other hand, Eq. (|5^) shows that the parametric suppression of the decay is largely 
overcome due to a numerical enhancement. The latter depends crucially on the particular 
values of the integers p,q. In particular, for our favorite choice p = IIW — 2Nf,q = 8, 
Eq.(^^ yields a factor ~ 10, while for p = Nc, g = 1 (as motivated by SUSY, see Sect.3) it 
is approximately two orders of magnitude larger, but still much smaller than the estimate 
of for the VVW potential. For a discussion of these results, see 
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To conclude this section, we would like to note that there also exist other domain 
walls interpolating between different local minima of the effective potential The 

surface tension a and decay rate T for these walls strongly depend on the vacuum states 
connected by the wall. 


6 Pseudo-goldstone bosons at different 9 angles 


In this section we address a few related questions. First, we discuss the calculation of the 
p' mass from the effective chiral Lagrangian and show that the main contribution 
to m^/ is given by the conformal anomaly. We also calculate the 6 dependence of the p' 
mass. Furthermore, it will be shown that for non-zero values of 6 , the pseudo-goldstone 
bosons cease to be the pure pseudoscalars, but in addition acquire scalar components. 

To study the properties of the pseudo-goldstone bosons, we parametrize the chiral 
matrix (^) in the form 


U 


Uq exp 




(55) 


®The Nc dependence displayed in Eq.(^) may look suspect as it apparently indicates that F = 
(P(g-i/Vc) Such a conclusion would be wrong, as in Eq.(^l|) we have neglected the second 

term in the effective potential in comparison to the first one. However, from the point of view of 
the Nc counting, the second term in Eq.(^^ is just the leading one. Therefore, it would be erroneous 
to extrapolate Eq.(Q) to the limit of very large Nc- One can easily check that in the limit Nc ^ oo 
with niq ^ 0 fixed, the lifetime of a metastable vacuum goes to infinity, in agreement with the picture of 
Witten Q for the case of pure gluodynamics. 
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where Uq solves the minimization equations for the effective potential (P^, and the helds 
7r“ , 7]' all have vanishing vacuum expectation values. A simple calculation of the matrix 
of second derivatives yields the following result for the mass matrix for an arbitrary value 
of 0: 
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{Mu cos (j)u + Md cos 0d) 

{Mu cos + Md cos 4)d + 4Ms cos (f)s) 


* IJ) I “K-p" ^ p ^ ^ t) + .5,““ * 
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{Mu cos (j)u - Md cos (pd) 


(56) 


mi3 = 2 


mi8 = 2. 



2 1 


NfUU 

~ 1 


3iV; UU 


{Mu cos (i)u - Md cos 0d) 


{Mu cos (pu + Md cos (pd - 2Ms cos 0*) , 


where <pi are solutions of the minimization equation stemming from Eq. (p7D . They depend 
on 9 as well as other parameters of the effective Lagrangian. If the 7i^ — rj — rj' mixing is 
neglected, mu coincides with the physical mass of the r]'. For 6 = 0 and the particular 
choice p = 3b,q = 8, we reproduce in this limit the relation given in |^: 



^ ^ rriii^i^i) +0{ml) . 

'f u,d,s 


( 57 ) 


(The choice p = N^, q = ^ would produce a numerically close result.) This mass relation 
for the p' appears reasonable phenomenologically. Note that, according to Eq. (pTl) , the 
strange quark contributes 30-40 % of the p' mass. This may lead us to expect that 
chiral corrections 0{m‘l) could be quite sizeable. We also note that in the formal limit 
Nc —>■ oo,mq —>■ 0 Eq. (|57D coincides with the relation obtained in |^. In this limit 
m^/ scales as Nf/N^, in agreement with Ref. In the different limit when W goes 

to inhnity at hxed non-zero rUg, the result is m^/ = 0{mqN^), as for ordinary pseudo- 
go Idstone bosons. As for the 6 dependence of the p' mass in the same limit, it is given by 
the third of Eq. (|5^ where the phases (pi implicitly depend on 9 through the minimization 
equation for the effective potential (Pif). 

The mass matrix (p^) can be used to study the mixing between pseudo-goldstone 
bosons (including also its 9 dependence). Let us consider the simplest case of the p — p' 
mixing which decouples from the in the isospin limit m^ = rrid = m, {dd) = (uu) = 
(qq). For the case 6* = 0, it is easy to verify that the mixing matrix 


1T7,q—.qi 


‘ -^(2m,(ss)+m(gg)) - m{qq)) 

4^''ms(ss) -m{qq)) ^ (§)^ 




{{as/TT)G^) 


K' 


(58) 
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coincide with an accuracy 0(mq) with the matrix given by Veneziano m 


m^_^, = 




2v^ l 
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mz 


+ \ml 
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X 

Nc 


( 59 ) 


with the only (but important) difference that the topological susceptibility ~ x in pure 
YM theory in the latter is substituted by the term proportional to the gluon condensate 
in real QCD in the former. For the particular values p = 3b = lliVc — 2Nf, g = 8, we may 
write down a QCD analog of the Witten-Veneziano formula: 

i^G^) = {ml + - 2ml) + . (60) 

This relation generalizes Eq. (|^ as it now includes the p — rj' mixing. 

Eqs. (^61) can also be used to study other problems related to the physics of the 
pseudo-goldstone bosons. In particular, we may find the mixing angles in the system 
at zero and non-zero angles 6. Instead of discussing these more phenomenological 
issues, we here would like to address another interesting aspect of the corresponding 
physics. Namely, we would like to show that the neutral pseudo-goldstone bosons in the 
^-vacuum cease to be the pure pseudoscalars, but instead become mixtures of the scalar 
and pseudoscalar states^ To show this, we note that the result found for the quark 
condensates in the 0-vacuum 


= ( 61 ) 

can be represented as a chiral rotation of the usual 0 = 0 vacuum: 

= {^u'^Rj)e=o • (62) 

Under such a rotation, the quark fields transform as 

«''« = (C'h)« 'T® . n'i = 'fu . (63) 

In the “rotated” basis, the spin content of the pseudo-goldstone bosons is the standard 
one. However, the relations (^) imply that it will generally have a different form in terms 
of the original unrotated fields. From Eqs.(|6l| j 6^ ) we obtain 

= (no sum over t) ^ (f/^),, = , (u+)^. = • (64) 

Let us consider e.g. the 7r° field in the 0-vacuum. In the ’’rotated” basis, it has the usual 
spin content. Using the correspondence (p5|J6^ , we obtain 

|7r°) ~ — d'ip^d!) = cos0„ \uip^u) — sin^^ \uu) — (u d) (65) 

Eq.(|55D illustrates the phenomenon announced in the beginning of this section: In the 
presence of a non-zero angle 0, the pseudo-goldstone bosons cease to be the pure pseu¬ 
doscalars, but in addition acquire scalar components. Although in reality 0 is extremely 
close to zero, this observation is not only of academic interest. The point is that in heavy 
ion collisions one can effectively create, in principle, an arbitrary value of 0 [Q. In these 
circumstances, the scalar admixture in the pseudo-goldstones would be quite large, and 
probably could play an important role in dynamics. 

® This fact was previously noted in the literature . 
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7 Further applications and speculations 


7.1 Axion potential from effective Lagrangian 

One of the interesting implications of the present effective Lagrangian approach concerns 
the possibility to construct a realistic axion potential |Q consistent with the known Ward 
identities of QCDQ. This may be achieved due to the one-to-one correspondence between 
the form of the axion potential V(a) and the vacuum energy E^acid) as a function of the 
fundamental QCD parameter 9. Indeed, the axion solution of the strong CP problem 
suggests (see e.g. [^) that 6 parameter in QCD is promoted to the dynamical axion 
held 6 —>• a{x)/fa, and the QCD vacuum energy E^aciO) becomes the axion potential 
V{a). Therefore, the problem of analysing V{a) amounts to the study of Eyac{0) in QCD 
without the axion, which is exactly the problem addressed above in the present paper. 

Both the local and global properties of the axion potential can be analysed with this 
approach. As for the former, we note that, as all dimensionful parameters in our effective 
Lagrangian are hxed in terms of the QCD quark and gluon condensate, the temperature 
dependence of the axion mass (and of the entire axion potential) can be related with that 
of the QCD condensates whose temperature dependence is understood (from lattice or 
model calculations). 

In particular, the axion mass, which is dehned as the quadratic coefficient in the 
expansion of the function E^ac{9) at small 9, is proportional to the chiral condensate: 
m\{T) ~ mq(0|TT|0)T//a- Therefore, m\{T) is known as long as (0|TT|0)r is known. 
This statement is exact up to the higher order corrections in rriq. We neglect these higher 
order corrections everywhere for T < T^. [T^. ~ 200MeD is the critical temperature), 
where the chiral condensate is nonzero and gives the most important contribution to m^. 
For the particular case Nf = 2 one expects a second order phase transition and, therefore, 
ml ~ ^(OlTTlO) ~ |Tc — T|^ for T near ~ 200 MeD. This is exactly where the axion 
mass does “turn on”. The critical exponent in this case [3 ~ 0.38, see e.g. recent reviews 


37| for a general discussions of the QCD phase transitions. 

The global (topological) structure of the axion potential appear to be rather compli¬ 
cated, in contrast to what could be expected according to simple model potentials such 
as V{a) = mla?/2 or V{a) ~ cos(a//a). In particular, it admits the appearance of ad¬ 
ditional local minima of an effective potential. Thus, the axion potential may become 
a multi-valued function, i.e. there would be two different values of the axion potential 
yi,2{9 = a/fa) for a hxed 9, which differ by the phase of the chiral held. For the VWW 
potential, this happens only at 6* ~ tt for a small isospin breaking, while for the ehec- 
tive potential (^^ with q ^ 1 metastable vacua exist for any 9, similarly to the case of 
a softly broken SUSY. Interpolating between two minima is the domain wall that was 
described in Sect. 5. We stress that it is not an axion domain wall, as the value of 9 does 
not change in this transition. Similar domain walls which separate vacua with diherent 
phases of the gluino condensate have been recently discussed p, |[ for the SUSY models. 
The appearance of such a domain wall implies an interesting dynamics which develops at 
temperatures below the chiral phase transition, which still has to be explored. 


Perhaps, one should note that some popular Ansdtze for the axion potential - such as V (a) = mfc? 
or F(a) ^ cos(a//a) - are at variance with the Ward identities of QCD. 
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7.2 Axially disoriented chiral condensate and axion search at 
RHIC 

Another interesting phenomenon amenable to an analysis within the effective Lagrangian 
framework is the possibility of prodnction of finite regions of 6 ^ 0 vacna in heavy ion 
collisions [Q, where the chiral helds are misaligned from the true vacuum in the axial 
17^(1) direction. This is somewhat analogous to the production of the disoriented chiral 
condensate (DCC) with a “wrong” isospin direction (see e.g. 


for a review) 


Let us briefly recall the reason as to why the DCC could be produced and observed 
in heavy ion collisions. The energy density of the DCC is determined by the mass term: 


= --Tr{MU + M+f/+) = -2m|(TT)| cos(0) 


( 66 ) 


where we put = rud = m for simplicity, and 0 stands for the misalignment angle. Thus, 
the energy difference between the misaligned state and true vacuum with 0 = 0 is small 
and proportional to rUg. Therefore, the probability to create a state with an arbitrary (j) at 
high temperature T ~ Tc is proportional to exp[—V{E^ — Eo)/T] and depends on 0 only 
very weakly, i.e. 0 is a quasi-flat direction. Just after the phase transition when (TT) 
becomes nonzero, the pion field begins to roll toward 0 = 0, and of course overshoots 
0 = 0. Thereafter, 0 oscillates. One should expect coherent oscillations of the n meson 
field which would correspond to a zero-momentum condensate of pions. Eventually, these 
classical oscillations produce the real tt mesons which hopefully can be observed at RHIC. 

We now wish to generalize this line of reasoning to the case when the chiral phases 
are misaligned in the Ha( 1) direction as well. For arbitrary phases 0* the energy of a 
misaligned state differs by a huge amount ~ E from the vacuum energy. Therefore, ap¬ 
parently there are no quasi-flat misaligned Ha( 1) directions among 0* coordinates, which 
would lead to long wavelength oscillations with production of a large size domain. How¬ 
ever, when the relevant combination (X)i 0* ~ from Eq.(^) is close by an amount 
~ 0{mq) to its vacuum value, a Boltzmann suppression due to the term ~ E is absent, 
and an arbitrary misaligned |6')-state can be formed. In this case for any 6 the difference 
in energy between the true l^) vacuum and a misaligned 16^)- state (when the 0* helds are 
not yet in their hnal positions 006*)) is proportional to niq and very small in close analogy 
to the DCC case. 

Once formed, such a domain with 6* 0 0 could serve as a source of axions, thus 
suggesting a new possible strategy for the axion search. Below we would like to sketch 
this idea, referring the interested reader to Ref. p3 for more details. 


It is well known that 0 is a world constant in the usual inhnite volume equilibrium 
formulation of a held theory. The superselection rule ensures that the only way to change 
6 under these conditions is to have an axion in the theory. However, due to the fact 
that we do not expect to create an equilibrium state with an inhnite correlation length 
in heavy ion collisions, the decay of a |6*)-state will also occur due to the Goldstone U 
helds with specihc CP-odd correlations^. Therefore, two mechanisms of the relaxation 

A similar phenomenon has been recently discussed in Ref.|^ where the possibility of spontaneous 
parity breaking in QCD at T ~ Tc was studied using the large Nc Di Vecchia-Veneziano-Witten effective 
chiral Lagrangian. 














of a |6')-state to the vacuum would compete: the axion one and the standard decay to 
the Goldstone bosons. In the large volume limit if a reasonably good equilibrium state 
with a large correlation length is created, the axion mechanism would win; otherwise, the 
Goldstone mechanism would win. In any case, the result of the decay of a |6*)-state would 
be very different depending on the presence or absence of the axion held in Nature^. 
Provided the axion production is strong enough, the axion could be detected by using 
their property of conversion into the photons in an external magnetic held [0. Thus, 
heavy ion collisions may provide us with a way to hnally catch the so far elusive axion. 


7.3 Early Universe during the QCD epoch 

As was discussed in Sect. 5, the ehective Lagrangian approach developed in this paper 
predicts the existence of the domain wall excitations in QGD at zero temperature. One 
may expect that these domain walls appear also for non-zero temperatures T < where 
Pc is the temperature of the chiral phase transition. If so, it would be very interesting 
to study this dynamics in the cosmological context. Here we only mention that the walls 
discussed above are harmless cosmologically as they decay in a proper time |^. On the 
other hand, as was noted in [Q , the dynamics of the decaying domain walls is an out-of- 
equilibrium process with 100% violation of CP invariance. This is because the phase of 
the chiral condensate in the metastable vacuum is nonzero and of order 1, that leads to 
violation of CP even if 0 = 0. (This is not at variance with the Vafa-Witten theorem 


which refers to the lowest energy state only.) It was speculated in that such effects 
could lead to a new mechanism for baryogenesis at the QGD scale. Indeed, it appears 
that all three famous Sakharov criteria could be satished in the decay of a metastable 
state discussed above: 

1. Such a metastable state is clearly out of thermal equilibrium; 

2. GP violation is unsuppressed and proportional to rriumdmsOeff, 9eff ~ 1. As is known, 
this is the most difficult part to satisfy in the scenario of baryogenesis at the electroweak 
scale within the standard model for GP violation; 

3. The third Sakharov criterion is violation of the baryon (B) number. Of course, the 
corresponding U(l) is an exact global symmetry of QGD. However, a “spontaneous” 
baryon number non-conservation could arise in this dynamics as a result of interactions 
of fermions with the domain wall. In this case, baryogenesis at the QGD scale is feasible. 
One possible scenario P| of such a “spontaneous” baryogenesis with zero net baryon 
asymmetry is a mechanism based on a charge separation, when the anti-baryon charge 
is concentrated on the surface of balls of the metastable vacuum produced in evolution 
of domain walls (B-shells). Rough estimates show that the observed ratio {ns — 
nB)ls r'U 10 ® can be easily reproduced in this scenario. Surprisingly, the energy density 
VLb associated with these B-shells can be close to unity. Therefore, they can be considered 
as candidates for dark matter. We would like to emphasize that each step in such a 
scenario for baryogenesis at the QGD scale could be, in principle, experimentally tested 
at RHIG. 


^^The possibility of production of axions in heavy ion collisions was independently discussed by Melissi- 
nos ffl. 
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Appendix: Fixing p/q by the “integrating in” 


The purpose of this appendix is to suggest a method which allows one to £x the number 


p/q that appears in the effective potential provided two plausible assumptions are 
made. One of them is insisting on the standard form of the fermion mass term in the 
effective potential, while the other one is the hypothesis of preserving the holomorphic 
properties when a heavy fermion is integrated in/out, see below. The approach developed 
below closely follows the method of where a similar problem was addressed for the 
effective Lagrangian for pure YM theory. The essence of this method is to consider 
QCD with light fermions as a low energy limit of a theory including in addition a heavy 
fermion, and to construct an effective Lagrangian for the latter theory starting from the 
effective Lagrangian ([I^). As will be shown below, a relation between the holomorphic 
and “topological” properties of two Lagrangians is non-trivial, and allows one to fix the 
crucial parameter ^ = q/{2p) entering Eg. ([T^) . 

The task of constructing such an effective Lagrangian for the theory with a heavy 
fermion is achieved by using the “integrating in” technique, developed in the context of 
SUSY theories in Ref, and reviewed by Intriligator and Seiberg The integrating 


in procedure can be viewed as a method of introducing an auxiliary held into the effective 
Lagrangian for QCD with light havors. Using the renormalization group properties of the 
QCD effective Lagrangian in the chiral limit rUg 0, the latter is extended to include 
the auxiliary held T, which will be later on identihed with the chiral combination QlQr 
of a heavy fermion. 

To conform with the notation and terminology of Ref.[^, we will call QCD with light 
quarks and the theory with a heavy fermion the d-theory (from “downstairs”) and the 
u-theory (from “upstairs”), respectively. The ehective potential of the d-theory is then 
Wd + Wj/ with (see Eq.(|T|)) 


Wd{h,U) = ]hLog 
4p 


h DetU 


S^qcd) 


^-ie 


- -Tr MU , 
2 


(A.l) 


(here c is a dimensionless numerical coefficient), and the summation over all branches 
of the logarithm in the partition function is implied. In this section Eq. (|A-l|) will be 
understood as representing a branch (section) of the multi-valued ehective potential, which 
corresponds to a lowest energy state for small 0 <C tt. As was shown in , this section 
corresponds to the principal branch of the rational function in the logarithm in Eq.( A.l ). 

We now want to relate |^, ^ the dimensional transmutation parameter Aqcd of the 


d-theory to the scale parameter Aqcd+i of the u-theory including a heavy quark of mass 
m ^ Aqcd, Aqcd+1- We assume both parameters to be dehned in the MS scheme, in 
which no threshold factors arise in corresponding matching conditions. The matching 
condition then follows from the standard one-loop relations 


Aqcd 

Aqcd+1 





bgcnYiMo) j 
SY 

&qcd+i5'^(Md) 


^QCD 


= h = —Nc--Nf 

3 3^ 


Ycd+1 — ^ Nc--{Nf+ 1) , (A.2) 
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and the requirement that the coupling constants of the d- and u- theories coincide at the 
decoupling scale Mq = m. We obtain 


^QCD — ^QCD+1 I A 2 
\ -/Vr 


m 


4/(3b) 


(A.3) 


QCD+l, 


As was explained in Ref.|]l7| Eq. (|A.3|) reflects the fact that, for fixed Aqcd+i, 

the scale parameter Aqcd characterizes the low energy theory surviving below the scale 
m, and thus depends on m. In this sense, the constant in the logarithm in Eq. ( |A.1| ) also 
depends on m 

/ \ 8p/{3bq) 

( T ) • (^-4) 

VAqcd+1 J 

Following Ref.[^, we now wish to consider (a particular branch of) the effective potential 


( A-l|) as the result of integrating out the auxiliary held T in the new effective potential 
W = Wu — mT which corresponds to the u-theory: Wd{h,m) = W{h,m, (T)), or 

Wd=[W^-mT](^T) , (A.5) 

where (T) is a solution of the classical equation of motion for the auxiliary held T : 

— m = 0 . (A.6) 


dT 


Let us note that, according to Eq. (A-5|) , Wd should depend holomorphically on (T). Our 
assumption is that this is only possible if an ehective potential W of the u-theory is itself 
holomorphic in the held T. Furthermore, one can see that Eqs. (|A-5|) ,( |Al^ ) actually dehne 
the potential Wd as the Legendre transform of Wu- Therefore we can hnd the unknown 
function Wu from the known potential Wd by the inverse Legendre transform: 


Wu = [Wd + mT](m) , 


where (m) solves the equation 


A 

dm 


{Wd + mT) = 0 . 


(A.7) 


(A.8) 


Eq. ([A.8[ ) can be considered as an equation of motion for the auxiliary “held” m. It is 
important to note that Eqs. ( |A.5| - |A.8|) imply that m should be treated as a com¬ 
plex parameter to preserve the holomorphic structure of Eq. (|A-l|) - When substituted 
in Eq. (|A-7|) , a solution (m) of Eq. (|A-8|) dehnes the potential Wu{h,T, (m)). When this 
function is found, the ehective potential W of the u-theory is dehned by the relation 


W{h, T, m) = Wu{h, T, (m)) - mT , 


(A.9) 


in accord with Eq. (|A.5|) . 

The solution of Eq.( |A.8| ) is easy to hnd using Eqs. 0,(0): 


2 h 
(m) = — — 
^ ' 3b T 


(A.10) 
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Thus, Eq.( |A.7D yields 


W^ = -\^h\og 

4 p 


VA4 \P/i 
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8p/{3bg) 
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QCD+l 
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DetU 


3b 


(ATI) 


Finally, Eq.(|A.9|) results in the effective potential of the u-theory 


W = 


1 g 
4 p 


h log 


(l-8/3fe)p/q 


T \ Det U 


. ‘^^QCD+l, 




:3—i6 


H —rh — mT 
3b 


(A.12) 

We expect this effective Lagrangian to describe QCD with light quarks and the additional 
heavy quark, corresponding to the held T, such that integrating out T brings us back to 
the effective Lagrangian ( AT ) for QCD. Indeed the equation of motion for the held T 
stemming from the ehective potential (|A.12|) reads 


m{T) = — h . 
' ' 3b 


(A.13) 


Inserting this classical vacuum expectation value (VEV) back to Eq. (|A.12|) (i.e. integrat¬ 
ing out the held T), we reproduce the ehective potential of the d-theory, Eq. ([A.l| ). Note 
that as Eq. (|A.13|) should preserve the Nc counting rule, we obtain {h) r'U Nc, 

(T) Nc- The Nc dependence of the VEV (T) is consistent with the identihcation 
(T) {QlQr) which will be suggested below. 

To identify the held T of the ehective theory with a corresponding operator of the 
fundamental theory, we note the following. As is seen from (|A.13|) , T has dimension 3, 
and thus should describe the VEV of an operator bilinear in the heavy quark helds. Fur¬ 
thermore, as long as m is ehectively considered as a complex parameter, this operator can 
only be QlQr or QrQl, in accord with the structure of the mass term in the underlying 
fundamental theory. Comparing Eq.( [A.13|) with the relation]^ between the VEV’s in the 
underlying theory: 

{mQLQn) = + i^Gg) , (A.14) 


247r 


and recalling the dehnition of the helds h, H, we conclude that 


(A.15) 


We thus see that the introduction of the heavy quark into the ehective theory hxes the 
parameter ^ which enters the ehective Lagrangian (|A.1|) for QCD. This has been obtained 
by the matching the local holomorphic properties of the d- and u-theories within the 
integrating in procedure. As was shown in |OT, the matching of the two theories can 


also be considered at the level of global quantization rules for the h , h helds which, for 
general values of p/g, are diherent for the theories described by Eqs. ( A.l ) and ( A.12 ). 


This matching agrees with the result ( |A.15|) , and provides a self-consistency check of the 


The equation ( A.14 ) follows from the operator product expansions {mQQ) = —(as/(127r)G^) 
0(l/m^) , {mQijsQ) = {as/{?>TT)GG) +0{l/rn^). 
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whole procedure of using the integrating in procedure to £x the parameter q/p = 2^ which 
enters the effective potential (^). Finally, we note that the correspondence 

mT , fhT . (A. 16) 

between the operators of the effective and underlying theories has the same meaning as 
Eqs.®, i.e. the classical held T describes the VEV of the chiral combination QlQr of 
the full (QCD +1) theory. 
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